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Extrapolating parton fragmentation to low Q 2 in e + -e collisions 

Thomas A. Trainor and David T. Kettler 
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We analyze the energy scale dependence of fragmentation functions from e + -e~ collisions using 
conventional momentum measures x p and £ p and rapidity y. We find that replotting fragmentation 
functions on a normalized rapidity variable results in a compact form precisely represented by the 
beta distribution, its two parameters varying slowly and simply with parton energy scale Q. The 
resulting parameterization enables extrapolation of fragmentation functions to low Q in order to 
describe fragment distributions at low transverse momentum pt in heavy ion collisions at RHIC. 
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I. INTRODUCTION 

QCD theory predicts that an abundance of soft gluons 
(minijets) should be produced in relativistic collisions at 
RHIC Copious gluon production may drive forma- 
tion of the colored medium in heavy ion collisions and 
global hydrodynamic phenomena (2, L|- However, the 
degree of equilibration of minijets in heavy ion collisions 
remains uncertain theoretically and experimentally. We 
should therefore search for and study remnants of low-Q 2 
(energy scale Q ~ 1 - 5 GeV) partons in single-particle 
and two-particle distributions of final-state hadrons. 

Measurements of two-particle correlations at RHIC 
with novel techniques have revealed substantial unequili- 
brated low-Q 2 _parton fragment structure in p-p and Au- 
Au collisions |illlllllal3llSIIllllllllli Frag- 
ment correlations have been measured in p-p collisions on 
transverse rapidity yt (defined below) and complemen- 
tary angular subspace (rj, <f>) with no jet hypothesis (no 
high-p t trigger particle), providing access to fragments 
from minimum-bias partons (no analysis constraint on 
parton momentum) dominated by minijets. Jet correla- 
tions have been observed in p-p collisions for hadron p t 
down to 0.35 GeV/c (parton Q ~ 1-2 GeV) HQ. Sim- 
ilar measurements in heavy ion collisions have revealed 
unexpected complexity. 

Initial measurements of two-particle angular correla- 
tions in p-p collisions at Fermilab (fixed target) and the 
CERN ISR on momentum subspace (77, (f>) (pseudorapid- 
ity and azimuth) [15J were described in terms of lon- 
gitudinal (string) fragmentation ^(|. Jets — correlated 
fragments from hard-scattered partons — were first ob- 
served at larger p t and y/s, establishing the nature of hard 
parton scattering and the validity of perturbative QCD 
(pQCD) [l^. A pioneering study of two-particle frag- 
mentation functions in LEP e + -e~ collisions is described 
in Q (c/. a related theoretical treatment in fl9jl. 

Minijet correlations in p-p and heavy ion collisions ob- 
served at RHIC represent QCD in a non-perturbative 
context: parton scattering, energy loss and fragmenta- 
tion at low Q 2 . We wish to connect those measurements 
to QCD theory via parton fragmentation measurements 
at larger Q 2 in elementary collisions. A context for two- 
particle fragment distributions at RHIC can be estab- 



lished by studying single-particle fragmentation functions 
(FFs) from p-p and e + -e~ collisions, the latter provid- 
ing especially precise access to the fragmentation process 
down to very low parton Q and hadron momentum. 

Low-Q 2 fragmentation is related to local parton-hadron 
duality (LPHD) which provides a correspondence be- 
tween pQCD parton predictions and hadronic observ- 
ables |20]. According to LPHD conversion of partons 
to hadrons occurs locally in configuration space, with al- 
most no distortion of the parton momentum distribution. 
The structural difference between a parton and a hadron 
should vanish for Q ~ 1 - 2 GeV where parton produc- 
tion is most abundant in RHIC collisions. LPHD is im- 
portant for low-Q 2 partons where 'fragmentation' may 
terminate with one or two partons (and hence hadrons). 
Minijet-related minimum-bias two-particle correlations 
studied in p-p and A- A collisions 0, 0, IE El EH 03 may 
therefore provide details of parton fragmentation at the 
energy scale where LPHD is most important. 

Modification of parton scattering and fragmentation 
in the QCD medium of heavy ion collisions may reveal 
properties of the medium itself. A theoretical study 
of in-medium modification of the single-particle FF in 
A-A collisions is reported in [2l|. The expectation is 
deformation of the in vacuo FF toward lower momen- 
tum, possibly corresponding to observed changes in the 
single-particle pt spectrum ratio Raa with collision cen- 
trality [22|, |23| . A related study of two-particle correla- 
tions in heavy ion collisions, especially the asymptotic 
approach of fragment distributions to thermal equilib- 
rium with increasing parton dissipation in the medium, 
is reported in |6(. Given the close connection between 
single-particle FFs and minijet-related two-particle frag- 
ment correlations, correlation measurements in heavy ion 
collisions may provide a more differential picture of prop- 
erties of the QCD medium and its influence on low-Q 2 
parton fragmentation. 

In this paper we establish a basis for extrapolation of 
measured e + -e~ FFs to small energy scales as prepara- 
tion for similar extrapolations in nuclear collisions. This 
is not a theoretical analysis based on pQCD methods. 
For a recent example of such an analysis which consis- 
tently describes FFs over a large x p (momentum fraction) 
range cf. [24| and the related discussion in Sec. IXIII Dl 
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The present paper describes a phenomenological analysis 
of FF data intended to provide the best possible extrap- 
olation down to small parton energies where pQCD as- 
sumptions such as collinearity and factorization become 
invalid. 

The paper is organized as follows: We first present a 
new method of analyzing fragmentation functions, with 
emphasis on rapidity y as a preferred kinematic variable 
for low-Q 2 fragmentation studies. We then consider the 
general properties of FFs in the context of the double-log 
approximation (DLA), angular ordering and color coher- 
ence. We compare measured FFs from e + -e~ collisions 
at three energies on several momentum variables and de- 
scribe a new form of approximate energy-scale invariance 
on normalized rapidity u. We demonstrate that FFs on u 
are precisely modeled by the beta distribution. We con- 
sider FFs for identified hadrons and identified partons. 
Based on fits to measured FFs and jet multiplicity data 
we develop a simple, precise parameterization of e + -e _ 
FFs valid over a broad energy range. Finally, we use our 
parameterization to study scaling violations and extrap- 
olation to low Q 2 . 



II. ANALYSIS METHOD 

The FF D(xe,Q 2 ) as used in this analysis is a single- 
particle density 2dn/dxE of hadron fragments on energy 
fraction xe = Ehadron/ 'Eparton produced by a pair of par- 
tons (dijet) with total energy Q (Q 2 = —q 2 is the negative 
invariant mass squared for the initial momentum trans- 
fer). Momentum fraction x p = Phadron/p P arton approxi- 
mates xe if particle momenta are measured. At large x 
the distribution shape reflects energy conservation during 
the parton splitting cascade. At small x the shape is de- 
termined by quantum coherence of gluon emissio n (g luon 
or color coherence and the hadron size scale) |25l l26|. 
The FF data used in this study are hadron distributions 
reported on momentum fraction x p , or £ p = ln(l/x p ). 
Distributions on x p emphasize pQCD aspects of parton 
fragmentation at large p (e.g., scaling violations). For 
non-pQCD effects, especially the role of gluon coherence, 
logarithmic variable £ p provides better visual access to 
the relevant small-Xp (large-£ p ) region. 

Dp (x, Q 2 or s) is the FF for parton type p and hadron 
type h at the energy scale denoted by Q 2 or s. The 
parton-flavor-inclusive distribution D h (x,s) is discussed 
in Sec. IVIII and the fragment-flavor-inclusive distribu- 
tion D p (x, s) is discussed in Sec. I Villi The total FF is 
D(x,s) = J^h D (x, s). The corresponding FF on £ is 
-0(£j s ) = xD(x,s), with Jacobian factor x. FFs sat- 
isfy relations J Q dxD(x,s) = 2n ( dijet fragment multi- 
plicity) and L dx x D(x, s) = 2 (dijet energy conserva- 
tion) |27| . To simplify notation we adopt the convention 
that symbol D represents any fragmentation function, 
with the specific form [Jacobian relation to D(x, s)] im- 
plied by the first argument. Plot axes are labeled with the 



corresponding dijet particle density D(x,s) — > 2dn/dx, 
-D(£, s) — > 2dn/d£, etc. to avoid confusion. 

This study focuses on low-Q 2 parton fragmentation. 
Since £ p and pQCD expansion parameter Y(Q) — 
ln(Q/A) (A represents a reference energy scale) are un- 
defined as p, Q — > we introduce rapidity y (well- 
behaved in that limit) as an alternative logarithmic mo- 
mentum/energy variable. The rapidity along axis z is 
y z (p; tuq) = ln[(E + p z )/mt], with transverse (to z) mass 
m 2 = tUq + p\. In frames comoving on z p — ► p tl 
E — > m t and y — > y t = \n{(m t + p t )/m a }. In a 
frame where p is the only non-zero momentum compo- 
nent y(p; mo) = hi[(E + p)/mo], with y — > ln(2p/rno) for 
p 3> mo and — > p/mo for p <C mo- mo may be a quark 
or hadron mass or energy scale A. 

Given the limiting cases for y we note that 
h^-y/i/mo) ~ y(y/s/2;m Q ) = y m ax, the kinematic limit 
for fragment rapidities. Similarly, Y(Q) = ln(Q/A) ~ 
y(Q/2;A) is a rapidity measure of the energy scale rel- 
ative to a reference scale. We observe for data a lower 
limit y m in which may depend on fragment species and 
collision system (e + -e~ vs p-p). For unidentified frag- 
ments we assign the pion mass mo — > m^ to all hadrons 
(but cf. Sec. IVIB . From data distributions on x p or 
£ p for parton energy scale Q or CM energy a/s we ex- 
tract fragment momenta p and calculate equivalent ra- 
pidities y (fragments) and y m ax (partons). Data dis- 
tributions on Xp or £ p are transformed to distributions 
on y using appropriate Jacobians. In |28l | ln(l/a;p) — > 

ln{ (E +p) parton /{E + p) hadron} = Vmax ~ V exactly. 

Most e + -e~ FFs plotted on normalized rapidity u = 
(y - y m in)/(ymax ~ y m in) ~ 1 - €p/Y have a partic- 
ularly simple form described by the beta distribution. 
The unit-normal beta distribution defined on u S [0, 1] 
is f3(u;p,q) — u p_1 (1 — u)^ 1 /B(p,q), with parameters 
p, q > and beta function B(p,q) = rfojff ■ Parame- 
ters p and q determine the shape of the distribution below 
and above the mode (most probable point) respectively. 
The mode is u* — p ~ 1 n , the mean is u = —7- and the 

p+q—2 ' p+q 

variance is a 2 u = { p +q y2 P { q p+q+1) ~ 4 (p+ 1 g +i) ( to 2% )- 



III. e + -e- FRAGMENTATION FUNCTIONS 

The double log approximation (DLA 25] ) provides 
a context for extrapolating the fragmentation process 
to low Q 2 . The fragment emission probability is ap- 
proximated by a uniform density on logarithmic space 
[log(0), log(p)], where 9 is the radiated-parton emission 
angle and p is its momentum. The distribution is 
sketched in Fig. ^ (left panel), where P is the leading- 
parton momentum and O ~ 1 is the jet angular accep- 
tance. The large solid triangle represents the kinematic 
boundary for the first radiated parton. The smaller tri- 
angles illustrate the self-similar nature of the splitting 
process (angular ordering |2|j), each radiated parton be- 
coming itself a DLA radiator. Alternatively, the DLA 
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may be expressed in terms of c?log(fc t ) dy, where k t is 
the transverse momentum component relative to the ra- 
diating parton momentum, and y is the radiated parton 
rapidity |30j . The flat DLA emission probability is ter- 
minated at some k t (g rey band in Fig. □ left panel) due 
to gluon coherence |25| . For sufficiently small k t the con- 
jugate transverse size of the virtual gluon overlaps the 
radiating parent parton and parton showering is termi- 
nated. 
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FIG. 1: Left panel: Schematic illustration of the double-log 
approximation (DLA) to parton fragmentation. Right panel: 
Self-similar variation with energy of the the corresponding 
fragment distribution on y. 

The general form of the FF and its evolution with 
Pparton corresponding to the DLA with angular ordering 
and gluon coherence (modified leading log approximation 
or MLLA |H3) is sketched in Fig. □(right panel). We ex- 
pect a monotonic increase with decreasing y below y max 
due to showering. The available phase space (above the 
band in the left panel) is however reduced with decreasing 
y by gluon coherence, causing the FF to turn over, with 
a maximum at y* (the mode) . The distribution then falls 
to zero at some y m in which may be nearly independent 
of y ma x- The FF is apparently self-similar at two levels: 
the internal cascade and its external boundary. Reducing 
the maximum opening angle (dash-dot line) or increas- 
ing the parton momentum P (dashed line) changes the 
boundary of the cascade (left panel) and correspondingly 
the FF (right panel). 



IV. FRAGMENT DISTRIBUTION ON x p AND £ p 

Single-particle FFs from e + -e~ collisions have been 
studied extensively (e.g., jH El HI S El E3)- 
The FF data plotted in Figs. [3- □ were obtained from 
collisions at three energy scales (CM energy Q = y/s 
= 14, 44 and 91.2 GeV) measured at PETRA |22| and 
LEP 26] for unidentified hadrons from unidentified par- 
tons (flavor- inclusive jets). Those distributions are fidu- 
cial for this study because of the exceptional data quality 
and fragment momentum coverage. We consider the data 
in several presentation schemes and then develop a pa- 
rameterized representation for extrapolation to low Q 2 . 
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FIG. 2: Left panel: 



tional momentum x p = 



: fragmentation functions on frac- 
Pfragmcnt /Pparton for three CM en- 
ergies. The dashed line is an exponential reference. Right 
panel: The same fragmentation functions on logarithmic vari- 
able £ p = ln(l/a;p). The vertical dotted lines mark equivalent 
points on the two variables. The solid curves are determined 
by the parameterization from this analysis. 



In Fig. |3 we plot fragment distributions on momen- 
tum fraction x p (left panel) and logarithmic equiva- 
lent £ p (right panel). Distributions on x p emphasize 
the large-Xp (small-£ p ) region where pQCD is expected 
to best describe data, where the naive parton model 
predicts 'scaling' or invariance of the parton distribu- 
tion on energy scale Q. The dashed reference line in 
the left panel illustrates the exponential model some- 
times used to characterize FFs on x p . The vertical dot- 
ted line corresponds to £ p = 1.5 in the right panel: 
only a small fraction of fragments fall above that point. 
The data exhibit systematic scaling violations (Q de- 
pendence) described by the DGLAP evolution equations 
(cf. Sec. HlHi. To study scaling violations FFs 
on x p arc parameterized by a model function such as 
D%(x, Q 2 ) = Nx a (1 - x f (1 + j/x), where the four pa- 
rameters depend on parton type p, hadron type h and 
energy scale Q |i(J, |41| . Distribution details in the small- 
x p region (e.g., below x p = 0.1) are minimized in this 
format (cf. Sec. lXLUDfl . 

We can also plot data on £ p (right panel) which em- 
phasizes the small-Sp (large-£ p ) region and better re- 
veals non-perturbative details of fragmentation. The dis- 
tribution is approximately gaussian, with mod e and 
r.m.s. width a^ p predicted by pQCD (cf. Sec. IXI|) . As 
noted, the fall-off at large £ c and maximum at £* result 
from gluon coherence |25|, |2(| . Measurement of the full 
fragment distribution above and below the mode is im- 
portant for a complete characterization of the fragmen- 
tation process. 

The solid curves in Figs. [21 -0 are obtained from beta 
distributions on normalized rapidity u determined by the 
systematic trends of beta parameters (p, q) plotted in 
Fig. [T2| (left panel) (cf. Sec. IX Cf) and transformed to 
each plotting space with appropriate Jacobians. Some 
approximation to 'scaling' or energy-scale independence 
is expected at large x p (small £ p ). Another form of seal- 
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ing at small x p (large £ p ) may be explored by plotting 
distributions on rapidity y. 



V. FRAGMENT DISTRIBUTION ON y 

Fragmentation functions plotted on £ p coincide at the 
kinematic limit £ p = corresponding to the parton mo- 
mentum. However, in Fig. (left panel) we observe that 
the FFs for three energies plotted on y have a common 
low-momentum limit y m in ~ 0.35 (vertical line, and cf. 
Sec. IIX All . That alignment is possible because y has the 
well-defined limiting value as momentum p — ► 0. Each 
data FF is terminated at the upper end by its kinematic 
limit y m ax = y(y/s/^', m o) (vertical lines) corresponding 
to £ p = in Fig. [5] The distribution maxima increase 
monotonically with collision energy. The FFs in the left 
panel illustrate the self-similarity sketched in Fig. ^ (right 
panel) and confirm an expectation for DLA scaling: frag- 
mentation at small y should be nearly independent of the 
leading parton momentum. 
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FIG. 4: Comparison of fragmentation functions for three 
energies on £ p and y m ax —y- Differences are noticible only for 
small y or p (large £ p ). 



Distributions on £ p extend in principle to oo, but the 
transformed beta distributions limit at p m in ~ m n /2 
or £ p ~ hi(Q/m n ), indicated by vertical lines in Fig. 0] 
(left panel). This comparison suggests that rapidity 



y(Q/2;mo) or difference y n 



y could replace £ p in 



FF studies. Rapidity y(Q/2; A) could also replace pQCD 
expansion parameter Y(Q) = ln(Q/A), remaining well- 
defined for Q — > while preserving established pQCD 
relations for larger Q. Figs. and 01 also suggest that 
rescaling the rapidity by y m ax — ymin might provide more 
differential access to FFs. 



FIG. 3: Fragmentation functions on rapidity y for e + -e 
collisions (left panel) and p-p collisions (right panel). 

In Fig. E3 (r ight panel) we plot FFs from p-p collisions 
at FNAL |33l |34j (the points are samples from the orig- 
inal data distributions used here to illustrate qualitative 
features). While the general features are similar to FFs 
for e + -e~ collisions the lower limit y m in is considerably 
larger for p-p collisions (note the dotted reference line 
common to the two panels). The larger y m in for p-p col- 
lisions (~ 1.5) may be due to the finite jet-cone ope ning 
ang le [34j and/or the presence of the underlying event [42j 
which must be distinguished from jet fragments. The 
gaussian curve labeled MB represents a minimum-bias 
fragment distribution (no selection is imposed on the 
parton momentum spectrum) derived from the event- 
multiplicity dependence of p-p pt spectra 01 which com- 
pares well with the systematics of FFs obtained from p-p 
jet reconstruction. 

In Fig.^we compare FFs on £ p = ln(Q/2p) (left panel) 
and y ma x — y ~ ha(Q/2p) (right panel). The distri- 
butions are equivalent below the upper half-maximum 
points (p 3> mo), above which distributions on y ma x — y 
drop rapidly toward well-defined limits at y max — y m in- 



VI. FRAGMENT DISTRIBUTION ON u 



Expectations of approximate energy scaling at large 
x p and a different form of scaling (gluon coherence) at 
small x p seem to require conflicting plotting strategies 
on £ p and y. However, both forms can be accommodated 
with normalized rapidity u = {y-y m in)/ (Umax ~ Vmin) € 
[0,1]. FFs from e + -e~ collisions can be factored as 
D{u,y max ) = 2n(y max ) g(u,y max ), with dijet multiplic- 
ity 2n(y max ) (cf. Fig. I13fl and unit-normal form factor 
g(u,y max )- In Fig. [S] we plot the three representative 
FFs transformed to l/n(y max ) dn/du = g[u, y m ax)- Mul- 
tiplicity 2n(y max ) can be obtained from fits to data, but 
also from the shape of g(u,y max ) (cf. Sec. IX B)l . 

We have determined that data form factor g(u, y ma x) 
is well-described by beta distribution (3(u;p,q) defined 
in Sec.[H] While there are substantial 'scaling violations' 
on Xp or £ p (cf. Sec. IXIl() . the normalized FF shapes on 
u in Fig. are nearly independent of Q 2 or y max over a 
substantial energy range. However, the remaining small 
variations with energy are significant, and well described 
by energy-dependent beta parameters (p, q) plotted in 
Fig. n~T1 and discussed in Sec. IX CI 
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FIG. 5: Fragmentation functions plotted on normalized ra- 
pidity u in linear (left panel) and semi-log (right panel) for- 
mats. The data distributions have been normalized by the 
corresponding di-jet multiplicity at each energy (lower solid 
curves in Fig. 1131 determined by parameters (p, q) . The data 
for three energies are plotted, but the curves for only 14 and 
91 GeV are plotted to provide visible separation. 



VII. IDENTIFIED HADRON FRAGMENTS 



To determine the effect of assigning the pion mass to 
unidentified fragments we used the following procedure. 
Data distributions on x p for three identified fragment 
species (the 91 GeV data in Fig. were transformed to 
normalized rapidity u with the proper mass assignments. 
Functions (3(u;p,q) were fitted to each species, trans- 
formed back to x p and plotted (dash-dot curves) with 
the data in Fig. (left panel). The model functions on 
x p were summed to represent the combination of uniden- 
tified hadrons and transformed to rapidity y assuming 
the pion mass, giving the solid curve in Fig. (right 
panel). The dotted curve was obtained by assigning the 
pion mass to all data, transforming to u, fitting the re- 
sulting distribution and then transforming back to y. We 
conclude from the results that misidcntifying kaons and 
protons as pions in unidentified hadrons shifts the FF 
peak mode at 91 GeV from the pion value ~ 0.38 to 
the inclusive hadron value ~ 0.41 in Fig. El (left panel). 
The dashed and dash-dot curves are fits to the individual 
fragment species with proper masses used to determine 
the rapidities. 



We define rapidities for unidentified hadrons by as- 
signing the pion mass mo — ► to several particle 
species. To assess the consequences we use identificd- 
particle FF data for two CM energies. In Fig. we show 
data g(u,y max ) and best-fit model /3(u;p,q) for identi- 
fied charged pions (left panel) and kaons (right 
panel) at 10 GeV and 91 GeV H^. Parton rapidity 
y-max is determined in each case with the identified hadron 
fragment mass. The distributions for identified protons 
(p,p) show similar behavior but with larger statistical er- 
rors. The pion FFs have widths similar to unidentified 
hadrons, but the peak modes are significantly lower (0.38 
vs 0.41 at 91 GeV). The kaon peak modes are comparable 
to those for unidentified hadrons but the peak width at 
higher energy is significantly larger. The kaon FF shape 
seems to converge on the pion distribution at lower en- 
ergy. The apparent blending of quark flavors could be 
related to the convergence of the gluon and quark FFs at 
lower energy in Fig. 1141 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



u u 

FIG. 6: Fragmentation functions for two CM energies and 
for pion (left panel) and kaon (right panel) fragments plotted 
on normalized rapidity u. 




FIG. 7: Left panel: Fragmentation function data for ^J~s = 
91 GeV and for three hadron species on fractional momen- 
tum x p and corresponding beta-distribution fits on normal- 
ized rapidity u transformed to x v (dash-dot curves). Right 
panel: Dashed and dash-dot curves for identified hadrons with 
proper mass assignments transformed from the left panel. The 
solid curve is a sum on x p of the beta-distribution fits from 
identified fragments in the left panel transformed to rapidity 
y with a pion mass assignment. The solid curve is close to the 
distribution obtained from inclusive hadrons (dotted curve) . 



From this exercise certain trends are notable: Pro- 
ton fragments have the largest momenta but the small- 
est rapidities. When transformed to normalized rapidity 
u the FFs for different fragment species are similar in 
shape (beta distribution) but exhibit small but significant 
mode variations with parton energy and hadron species 
(c/. Fig. [5] and the discussion following for a summary 
of flavor dependence). Unit-normal data distributions 
g(u, Umax) for all light hadron species are well-described 
by model /3(u;p, q), establishing applicability of the beta 
distribution to FFs for identified light meson and baryon 
fragments as well as to inclusive hadrons. 
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VIII. IDENTIFIED PARTONS 

We now consider the role of parton identity in FF 
systematics. Normalized data distributions on u are 
shown in Fig. [S] for inclusive hadrons from udsc-quark 
jets (upper- left) , gluon jets (upper right) and b-quark jets 
(lower-left) for several parton energies in each case [4|| 
l47j . The measured FFs for light quarks and gluons 
are well described by model f3(u;p, q), shape parameters 
(p, q) depending on parton species and energy scale. Di- 
jet multiplicities are obtained as the best-fit coefficients 
of the unit-normal beta distribution. As expected, there 
is a substantial difference between quark and gluon FFs 
at larger jet energies, and a strong energy dependence 
of gluon jet shapes for smaller jet energies evident in 
the upper-right panel (the two solid curves correspond 
to - 5 and 40 GeV gluons) (cf. Scc.ETJ. The b-quark 
data in the lower-left panel are not well described by 
the beta distribution. The best-fit beta distributions for 
= 91.2 GeV udsc quark and Q = 80.2 GeV (equiva- 
lent dijet energy) gluon jets (data with the best statistics) 
are repeated as the dashed (/3 q ) and dash-dot (/3 g ) curves 
respectively in all three panels to provide references. 




1 



i • 
J 1 " 

5 



5 ,„- 2 



10 
10 
10 







- 




b quarks ^qi\ 
















-H 






I A 


8.0-18.0 GeV 




□ 


18.0-38.0 GeV 




■ 


38.0-60.0 GeV 




O 


91.2 GeV 




1 ' 


183-209 GeV 












91.2 GeV udsc and b quarks 
80.2 GeV gluons p-°- 




0.25 0.5 0.75 



FIG. 8: Unit-normal fragment distributions on normalized 
rapidity u for udsc quarks (upper left), b quarks (lower left) 
and gluons (upper right) fragmenting to inclusive hadrons. 
Lower-right panel: Normalized logarithmic variable v(f3) (see 
text) measures the FF shape for b-quark jets (solid) relative to 
those for gluon jets (dash-dot) and udsc-quark jets (dashed) 
as references. 

In Fig. [S] (lower-right panel) we compare FFs from dif- 
ferent parton types in a more differential format. As 
noted, /3 q (dashed curves) for light fragments in udsc 
jets and /? g (dash-dot curves) for gluon jets are approx- 
imate limiting cases for all f3(u;p,q). We therefore de- 



fine v max = ln(/3 q + /3 g ), v min = -ln(l//3 q + l//3 g ) and 
normalized variable v([3) = (ln/3 - v mm )/ (v max - v min ), 
with v(f3 q ) + v(f3 g ) = 1. We plot v(f3 q ) (dashed), v(f3 g ) 
(dash-dot) and w(/3b) (solid) in the lower-right panel with 
the corresponding data for y/s = 91.2 GeV quarks and 
Q = 80.2 GeV gluons (also cf. data and solid curves in 
Fig. EH). 

The light-fragment distribution from b quarks (solid 
dots) coincides with v(f3 q ) (and open circles) for u < 0.7, 
but diverges sharply from the quark-jet trend above that 
point and descends towards v((3 g ) (and open triangles) 
for u > 0.7. The b-quark fragment data were reduced 
by 10% to coincide with the quark-jet curve below u ~ 
0.7. The initial normalization is represented by the beta- 
distribution fit v(Ph) (solid curve) with mode near 0.5. 
With this more differential format we confirm that fa- 
quark light-hadron fragments are not well described by 
a beta distribution. The exceptional softness of the fa- 
quark FF (for unidentified fragments) was anticipated 
theoretically (and cf. Fig.EJ. 




FIG. 9: Distributions on u for several quark/meson flavor 
combinations, showing evolution of the g(u, ymax) shape with 
quark/meson mass. The c — > D data are from |49j, and the b 
— * B data are from [5(J. The energies are dijet energies. 

In Fig. |5| we summarize FF data and models for sev- 
eral fragment and parton types. The pion, kaon and pro- 
ton FFs are beta-distribution fits to 91 GeV identified- 
fragment data (pion and kaon data are shown in Fig.|SJ|. 
The gluon FF is the beta distribution defined at 80 GeV 
by (p, q) systematics in Fig. 1141 (cf. comparison with FF 
data in Fig. ED ~ right panel) . The solid dots are b — > 
B data from [50| compared to a best-fit beta distribution 
(dash-dot curve) and theory (solid curve). Low-statistics 
c — * D data from [4^ are summarized by a best-fit beta 
distribution (dash-dot curve) and theory (solid curve). 

The two solid curves on the right of Fig. are from 
a theoretical treatment of heavy-quark fragmentation in 
which the FF for Q — > H(Qq) + q is approximated 
by D^(x p ) oc 1/{1 - l/x p - cq/(1 - x p )} 2 /x p , with 
eg cx 1/toq [HJ. The agreement of D^(u) (right-most 
solid curve) with b-quark data (solid points) [5(| us- 
ing 6q — 1.16/t7?| = 0.055 and m b ~ 4.6 GeV/c 2 is 
good. The dash-dot curve is the best-fit beta distribu- 
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tion with (p, q) = (23, 3) which does not describe the 
b — > i? data. That failure may be related to the ex- 
ceptional behavior of b — > light hadrons discussed in 
connection with Fig. [S] (lower panels). The solid curve 
for c — > D is D®(u) from the heavy-quark theory treat- 



ment, with cq = 0.57 /m 2 . = 0.29 and 



1.4 GeV/c 2 



The associated dash-dot curve is a beta distribution with 
(p, = (7.0, 2.8) which best describes the data from |49| . 
Both curves are consistent with the data, but the data 
errors are large below the FF peak mode. 

FF modes increase monotonically with increasing me- 
son and parton mass. However, the proton FF mode 
for udsc jets is lower than the inclusive hadron mode 
for gluon jets and the FF is significantly broader. The 
kaon FF shows the effect of the heavier s-quark mass, 
consistent with the trend for charm and bottom quarks 
(however, see the next paragraph). The FF mass de- 
pendence on normalized rapidity u is subtle compared to 
the kinematic dependence on meson and parton masses 
encountered on p u £ p or y. 

To summarize flavor dependence, the beta distribution 
describes the FF data for identified light quarks and glu- 
ons fragmenting to identified light mesons or baryons, 
providing a compact representation of the flavor depen- 
dence of fragmentation. The quality of the description is 
not good for heavy quarks fragmenting to light or heavy 
mesons. However, the region near u = 1 can be compared 
with non-perturbative trends for light-quark fragmenta- 
tion extrapolated to small Q 2 in Fig. [21 (right panel), 
where the fragmentation 'cascade' is a single splitting or 
no splitting (parton — > hadron). 



IX. FITTING f3(u;p,q) TO DATA 

We now fit the beta distribution to a sample of mea- 
sured FFs falling in three groups: 1) the five fiducial 
FFs for unidentified fragments from flavor-inclusive par- 
tons distinguished by nearly complete coverage of the 
kinematically-allowed fragment momenta [26t |3^ an d a 
selection of data for 2) identified fragments and 3) iden- 
tified partons to explore the role of hadron and parton 
species in fragmentation. Data in the form D(u,y max ) 
are fitted with model function 2n(y max ) [3(u;p,q) (cf. 
Sec. [nj, minimizing x 2 while freely varying parameters 
2n, p, and q, with u = (y - y m in)/{Vmax ~ Vmin) and 
Vmin constrained to specified values based on systemat- 
ics studies. 



A. Inclusive fragments from inclusive partons 

We first fit FF data for inclusive hadrons and par- 
tons. Table [I] contains the best-fit parameters for the 
five fiducial FFs (OPAL and TASSO [12 data) with 
Vmin = 0.35 (p = 0.05 GeV/c). The model functions 
with starred energies are compared to data in Figs. [21 - 
[SJ As noted, the FF shape is nearly independent of y/s, 



but there is a significant trend for q to increase and p 
to decrease with increasing energy scale, shifting the FF 
mode to smaller u. The fitted multiplicities agree with 
the q-q multiplicity curves and data in Fig. 1131 



TABLE I: Unidentified fragments from unidentified partons: 
Beta-distribution parameters from \ 2 fits to fragmentation 
functions for five energies. The FFs for starred energies are 
plotted in Figs. 0-0 



v^(GeV) 


2n 


P 


q 


X> 


14* 


8.8±0.10 


2.95i0.08 


3.52i0.07 


16/18 


22 


10.7L0.15 


2.91i0.06 


3.52i0.08 


25/20 


35 


13.4i0.05 


2.84i0.02 


3.50i0.02 


148/22 


44* 


14.6i0.10 


2.89 i0.03 


3.52i0.04 


49/22 


91.2* 


20.4i0.05 


2.84i0.01 


3.67i0.01 


86/51 




0.2 0.4 0.6 0.8 1 
u 



2 0.25 0.3 0.35 0.4 0.45 0.5 



FIG. 10: A study of the systematic effect on fiducial FFs and 
their beta parameters (p, q) of variations with y m in ■ The low 
ends of the FFs are most sensitive, and therefore p and FFs 
for low y ma x are affected most. The trends in the lower-right 
panel may be compared with the (p, q) entries in Table [I] 

Fig. shows the systematic dependence of beta pa- 
rameters (p, q) on the choice of y m i n ■ The first three 
panels illustrate the variation of fits and data on u with 
Vmin for energies 14, 44 and 91 GeV. x 2 variations are 
small over the interval shown. Variation with y m i n is 
greater for smaller u and smaller y ma x ■ Those trends are 
reflected in the summary of (p, q) variations in the lower- 
right panel: the p variation is greater, and more so for 
lower energy. The lines have slopes 0.8 (solid) and 1.6 
(dashed). We set y mi „ = 0.35 for all inclusive fits and 
discuss the related systematic uncertainties in Sec. IX CI 
The shift of OPAL (rescaled) q data in the lower-right 
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panel (closed squares to open circles) results from in- 
creasing all particle momenta by 6.5% to test the effect 
of uncertainty in the momentum calibration. The p data 
are much less affected. 



B. Identified fragments from inclusive partons 

We next explore the role of hadron identity in fragmen- 
tation, with mo — > m-hadron assigned for both fragment 
and parton rapidities. Fits to identified hadron fragments 
from flavor-inclusive partons at y/s = 10 GeV [44j and 
91 GeV plotted in Fig. Eland Fig. □(left panel) are 
presented in Table ITU Trends for fragmentation to light 
hadrons are summarized in Fig. [5] For pions, parame- 
ter p is smaller and q larger than for inclusive hadrons 
in Table UJ shifting the peak mode to smaller u as noted 
previously and understood as an effect of misidentifying 
kaons and protons as pions in the inclusive hadron frag- 
ment mixture. The pion fit x 2 is large; however the fit 
residuals are generally point-to-point random and sub- 
stantially larger than the stated errors, especially toward 
the ends of the distribution. 

TABLE II: Identified fragments from unidentified partons: 
Beta-distribution parameters from \ 2 fits to fragmentation 
functions for pions, kaons and protons at y/s = 8-18 and 
91.2 GeV 



FID 


2n 


P 


q 


ymin 


X> 


y/s = 8 - 18 GeV 




5.63±0.02 


2.92±0.03 


3.96±0.05 


0.35 


89/49 


K± 


0.88±0.025 


3.15±0.12 


3.84±0.19 


0.10 


33/39 


P, P 


0.18±0.02 


2.60±0.40 


4.30±1.00 


0.05 


21/24 


yfi = 91.2 GeV 


incl. 


20.4±0.05 


2.84±0.01 


3.67±0.01 


0.35 


86/51 




17.36±0.03 


2.66±0.01 


3.77±0.01 


0.35 


483/36 


K± 


2.39±0.03 


2.58±0.03 


2.99±0.01 


0.10 


10/26 


P, P 


1.10±0.02 


2.36±0.04 


3.58±0.07 


0.05 


17/23 



The 10 GeV kaon peak is similar to the pion peak. 
However, the 91 GeV kaon peak is much wider (p + q 
is reduced) and the mode is shifted substantially to the 
right (q — p is reduced) relative to the pion peak, consis- 
tent with the quark mass-dependence trend in Fig. [5] 
The proton peak mode is shifted further to the left, 
beyond the pion and gluon peaks as shown in Fig. [5J 
mainly by reduction of p. We note in passing that 

, |" toq+50 McV/c 2 " 



C. Inclusive fragments from identified partons 

Finally, we consider data for unidentified hadron frag- 
ments from identified partons for two parton classes 
shown in Fig. [5] udsc quarks (in combination) and glu- 
ons. The fit results are shown in Fig. 1111 and Table ITTT1 



The 'inclusive' table entry (first row) repeats the 91 GeV 
results from unidentified hadrons in Table [TJ for reference. 
Parameters for the free x 2 fit to the udsc FF data in the 
second row of the table reflect a width similar to the in- 
clusive data (q + p is similar), but the mode is shifted to 
slightly larger u (q — p is smaller). Details of the fitting 
procedure are shown in Fig. 1111 (left panel). The points 
are substantially larger than the reported errors. The free 
fit (solid curve) is strongly influenced by the single point 
at ii ~ 0.13. The \ 2 ls large , and the fit function misses 
the data near the peak. As for the pion fragment FF 
data the udsc data errors appear to be underestimated. 

TABLE III: Unidentified fragments from identified partons: 
Beta-distribution parameters from \ 2 fits to fragmentation 
functions for udsc quarks with yfs = 91.2 GeV and gluons 
with Q = 80.2 GeV. 



PID 


2n 


P 


q 


ymin 


x7" 


incl. 


20.4±0.05 


2.84±0.01 


3.67±0.01 


0.35 


86/51 


udsc 


18.36±0.04 


2.99±0.015 


3.55±0.015 


0.35 


209/19 


udsc 


param. (p, q) 


2.85±0.05 


3.58±0.05 


0.35 


550/19 


gluon 


27.2±0.4 


3.50±0.10 


5.10±0.15 


0.35 


5.7/22 


gluon 


param. (p, q) 


3.43±0.10 


5.30±0.10 


0.35 


14.4/22 



The dashed curve is constrained by the (jp, q) energy 
systematics in Fig. (left panel) consistent with fits to 
the fiducial inclusive FFs (second udsc row of the table). 
The peak of the data FF is better described, but the in- 
creased deviation from the small-u point greatly increases 
the x 2 - We expect the udsc FF to differ slightly from 
the flavor-inclusive FF (small shift to the right) due to 
the absence of gluon fragments, and the free fit is consis- 
tent with that expectation. We expect the udsc/inclusive 
multiplicity ratio to be 0.94 54] but observe 18.36/20.4 
= 0.90, suggesting that the udsc FF height is underesti- 
mated by the free fit. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



u u 

FIG. 11: Details of beta distribution fits to fragmentation 
functions for identified partons fragmenting to unidentified 
hadrons. The solid curves are free fits to data. The dashed 
curves are determined by the (p,q) systematics in Fig, 1141 (left 
panel) . 

In Fig. ^] (right panel) we show a free fit to FF data 
from gluon jets (solid curve) producing the fit parame- 
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ters in the first gluon row of Table ITTT1 which are plotted 
as open squares in Fig. ^] (left panel) and provide con- 
straints on the gluon (p, q) energy systematics discussed 
in Sec. IX Cl The free fit has an unusually small x 2 ; the 
data errors above the mode seem large compared to the 
residuals there. The dashed curve is a 'fit' with parame- 
ters constrained to the 80 GeV (p g , q g ) systematic values 
from Fig. ^] (left panel) and reported with its % 2 value 
in the second gluon row of Table IIIII 



X. ENERGY SCALE DEPENDENCE 

We can combine fits to fiducial FF data and dijet multi- 
plicity data to determine the energy dependence of (p, q) 
for quark and gluon jets over a broad energy range. Fits 
to data g{u,y max ) with model f3(u;p,q) determine spe- 
cific values (p, q) which constrain parameterized curves 
(p(y m ax),q(y m ax))- Fits to 2n(y max ) data via the (x E ) 
integral of f3(u; p, q) also constrain the parameterizations, 
especially important in energy intervals where there are 
no FF data available. The resulting (p, q) energy trends 
efficiently represent e + -e~ FFs over a broad energy range 
and provide a basis for extrapolating FFs to low Q 2 . 



A. Energy conservation sum rule 

The total FF D(x E ,s) = J2h D ' l ( x E,s) (sum over 
all hadron species) integrates to total dijet multiplic- 
ity ^mj/^^fi D(x E , s) = 2n to t{s) and satisfies the en- 
ergy sum rule (ESR) f 2mo /^dx E x E D(x e , s) = 2 [27|. 
The ratio of the integrals defines mean energy frac- 
tion (x E ) = l/n tot (s). Switching to (u,y,y max ), since 
D(u,y max ) = 2n tot (y ma x) g(u,y ma x) and g(u,y max ) ps 
f3(u;p,q) we have (x E ) w /„ dux E (u,y max ) (3(u;p,q), 
with x E (u,y max ) = cosh[y(ii)j/cosh(j/ max ) and y(u) = 
uy m ax + (1 — u)y rn i n . Those relations connecting 
/3(u;p, q) to n to t(ymax) are used below to obtain the en- 
ergy dependence of (p, q) from multiplicity data. 

Given several hadron species h with FFs D h (x E ) 
and dijet multiplicities 2nu we expect ESR = 
~^2 h 2rih (x E )h = 2, provided all species are integrated. 
However, if only charged hadrons are detected we expect 
ESR ~ 2/3 x 2 ~ 1.33. We can test the charged-fraction 
ESR using the fits to charged pion, kaon and proton 
data at 91 GeV from Table HJJ The ESR for inclusive 
FFs can be tested with the fits from Table [I] In general, 
if / is the ESR fraction for detected particles (/ ~ 2/3 
for the charged-hadron fraction) we expect the relation 
(x E ) = f/n(y max ) between monojet (charged-particle) 
multiplicity and energy fraction, which we use below to 
relate multiplicities to beta parameters (p, q) . 

Beta distribution fits to FFs g h (y or u, y max ) for iden- 
tified pions, kaons and protons at 91.2 GeV from Ta- 
ble [n] and the FF for inclusive hadrons are plotted on 
rapidity y and normalized rapidity u in Fig. 1121 The pa- 




FIG. 12: Unit-normal model functions @(y or u;p,q) from 
fits to FFs for identified fragments (pions, kaons and protons) 
and inclusive fragments, all from inclusive partons at y/s = 
91.2 GeV, plotted on rapidity y and normalized rapidity u. 
LONGER LINES 

rameters for those curves are used to obtain (x E )h for 
each hadron species using the correct hadron mass and 
(x E )inci for the inclusive distribution assigning the pion 
mass to all hadrons. We use the 2n fit values in the ta- 
bles to obtain J2h 2n h {x E ) h = 1.25 ± 0.03 (/ = 0.62) 
and 2n inc i (x E ) inc i = 1.18 ± 0.05 (/ = 0.59) for identi- 
fied and inclusive charged fragments. If (x E )i nc i is calcu- 
lated with the weighted-mean mass 0.2 GeV (weighted by 
the hadron multiplicities in Table [H| we obtain ESR — 
1.4 ± 0.05. The same procedure applied to the fits to 
lower-energy data from Table ITU gives ESR ~ 1.1. The 
exact energy scale for the lower-energy sum rule is not 
clear because of the scale range, but the result is roughly 
consistent with expectations. For the inclusive analysis 
below assuming the pion mass we use ESR factor 1.18. 

B. Dijet multiplicities from /3(u;p,q) shapes 

Dijet multiplicity 2n can be obtained directly by in- 
tegrating measured and extrapolated FF data, as in Ta- 
bleslTI- ITrTl However, as we have just shown there is a cor- 
respondence between 2n(y max ) and the shape of data FF 
g(u;y max ) or fitted model function /3(u;p,q) determined 
by parameters \p(y ma ,x),q(y m ax)]- We have obtained for 
inclusive charged fragments with pion mass assignment 
the relation 2n(y max ) = 1.18/ f Q dux E (u,y max )[3(u;p,q) 
at 91.2 GeV which we now use to relate energy trends 
of FF shape parameters (p, q) to fragment multiplicities. 
Measured multiplicity trends on parton energy thereby 
provide constraints on the energy dependence of FF pa- 
rameters (p,q), even in energy intervals where there are 
no measurements of FFs. 

Fig-ED s ho ws dijet multiplicities 2n for g-g and q-q par- 
ton pairs. Precise multiplicity data for quark jets from 
two-jet events have been available for some time. New 
methods have produced similarly precise gluon-jet multi- 
plicities from three-jet e + -e~ events. Data for gluon jets 
were obtained from CDF (closed triangles) j35j, CLEO 
(open triangles) |52(, OPAL 'jet-boost' algorithm (open 
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circles) [36| and OPAL inclusive (star) [37] . Data for 
quark jets were obtained from a compilation (Table 6 
in 1^3) and multiplied by factor 0.94 (the fraction of udsc 
jets in a flavor- inclusive sample [54]) to compare with 
the gluon jet multiplicities. The large points labeled ir 
and K are multiplicities from fits to identified fragment 
data HHHHIl in Table |Il| plotted with the indicated 
multipliers. The hatched regions represent the domain of 
low-Q 2 partons which motivates this extrapolation study. 
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FIG. 13: Dijet charged-particle multiplicity vs energy scale 
Q (dijet energy) plotted in a conventional format (left panel) 
and vs parton rapidity assuming the pion mass (right panel) . 
The solid curves are quark and gluon dijet multiplicities 2n q 
and 2n g obtained from the (p, q) parameterizations in Fig. 1141 
(left panel). The dash-dot curve in the left panel is from a 
3NLO pQCD expression. The udsc quark-jet multiplicities 
for unidentified hadrons (solid dots) are taken from a sur- 
vey in 53]. The dotted curves in the right panel illustrate 



quadratic trends A (y„ 



ymin) (see text). 



The solid curves in Fig. 1131 are multiplicities derived 
from the (p, q) energy trends using the relations defined 
above. The (p, q) parameterizations are adjusted to fit 
the multiplicity data but constrained by (p, q) values 
from fits to fiducial FFs. The resulting (p, q) energy 
dependence is described in the next subsection. Be- 
cause 2n oc 1/(xe) and (i^) is monotonic with mean 
u = p/(p + q), multiplicities are mainly determined by ra- 
tio q/p or difference q — p (i.e., the mode or mean of the 
beta distribution), and only weakly dependent on sum 
q + p (the width). A unique description of (p, q) over 
a broad energy range requires fits to multiplicity trends 
supplemented by the fits to fiducial FFs described in the 
previous section. 

Quark-jet multiplicities are described in the 
MLLA by 3NLO expression n q (Y) = if/2.25 • 

Y- aiC ' 2 exp{2CVY + ad{Y)}, with Y = ln(Vs/A), 
C = y/4n c /b and b = (lln c - 2n/)/3 56]. We used 
at = 0.3 from [E3| and K = 0.13 and A = 0.15 GeV 
from Tabic 1 and the functional form of 8 g (Y) from Fig. 
3 in [5(j for S(Y) (all for n/ = 5). We set the coefficient 
of S(Y) to a = 1.8 to obtain the best agreement with 
quark-jet data, shown by the dash-dot curve in Fig. 1131 
(left panel) just visible relative to our parameterization 
(solid curve). 



Variation of dijet multiplicities in the form A (y m ax — 
Urnm) 2 would be expected for the self-similar scaling il- 
lustrated in Fig. ^with fixed FF mode it*. Quadratic 
trends for quark and gluon jets are illustrated by the 
dotted curves in Fig. ^] (right panel), with A = 0.5 
for quarks and 1.45 • 0.5 for gluons. Deviations from 
the quadratic trend for quark-jet multiplicities in Fig. 1131 
correspond to the linear variation of (p q ,q q ) with y max 
above y ma x = 4.5 in Fig. ^] (left panel) which shift u* to 
smaller values, as illustrated in the right panel of that 
figure. Gluon-jet multiplicities deviate more dramati- 
cally from the quadratic trend at lower energies, mov- 
ing from the quark-jet curve to a gluon-jet trend about 
50% larger within the energy interval y m ax = 3.5 — 5 
(Q = 5 — 20 GeV) as the quark-gluon color charge differ- 
ence emerges. Above 20 GeV the gluon-jet multiplicities 
reflect the smaller linear variation of (p g , q g ) with y ma x 
in that energy interval. 



C. Energy dependence of (3(u\p,q) parameters 

Fig- E] (left panel) shows the (p, q) energy dependence 
which produces the quark and gluon jet multiplicities 
(solid curves) in Fig. ^| and the solid curves compared 
to fiducial FFs in Figs. [21 - EI Those curves summarize 
the energy dependence of udsc and gluon fragmentation 
to unidentified hadrons in e + -e~ collisions. We assume 
that the shapes of inclusive (dominated by light quarks) 
and udsc FFs are approximately the same, as in Fig. 6 
of |4y| . Inclusive and udsc jet FFs at 91.2 GeV are com- 
pared in Table II I II and Fig. ^2 (left panel) and found 
to be similar. The vertical dotted lines mark the limits 
of multiplicity measurements, while the vertical dash-dot 
lines mark the limits of measured FFs used in this anal- 
ysis. The upper ten solid points represent the fiducial 
FFs. The open squares represent the single gluon FF in 
Fig. 1111 fright panel) which constrains (p g ,g g ). 

The (p, q) curves in Fig. [21 (left panel) are described 

by 



Pq 

Ps 



2.90 ± 0.05 - (0.05 ± 0.01) (y max - 5.3) (1) 
3.50 ± 0.05 + (0.05 ± 0.01) {y max - 5.3) 
(0.8 ± 0.2) (y max - 4.5)[tanh(y ma;c - 3.5) - l]/2 
p q + [(0.07 ± 0.02) (y max - 5.3) + 0.55 ± 0.05] x 



{tanh[(2.5 ± 0.5) (y 

ma 

q q + [(0.07 ±0.02) (y„ 
{tanh[(2.5 ± 0.5) {y ma 



-4.1±0.1)] + l}/2 
- 5.3) + 1.70 ±0.07] x 
-4.1±0.1)] + l}/2 



The p q expression in Eq. JQ) is determined only by linear 
interpolation and extrapolation of fits to the hadron- and 
parton- inclusive fiducial FFs. Given that definition of p q 
the curve for q q is then defined only by the fit to the udsc 
quark-jet multiplicity data in Fig. 1131 The agreement in 
Fig. ll4l between q q determined by fitting light-quark mul- 
tiplicities 2n q (upper solid curve) and by fitting individ- 
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ual inclusive-parton FFs (upper solid points) indicates 
the consistency of the two methods. The expression for 
p g is guided by the fit to a single gluon FF denoted by 
the lower open square point in Fig. El (left panel) , but is 
also influenced by its impact on Fig. El The expression 
for q g is then determined relative to p g mainly by gluon 
multiplicities 2n g , but is also influenced by its impact on 
Fig- El See the discussion of that figure for more details. 




FIG. 14: Left panel: Beta-distribution parameters (p q ,q q ) 
and (q s ,p s ) respectively for light-quark (solid) and gluon 
(dashed) jets and corresponding gluon-to-quark-jet multiplic- 
ity ratio r vs parton rapidity y m ax ■ Right panel: Unit- normal 
FFs (beta distributions) obtained from parameters in the left 
panel plotted on normalized rapidity u for quark (solid) and 
gluon (dashed) jets and for nine equal-spaced values of parton 
rapidity y m ax illustrating peak shape evolution with energy 
scale. 

The error bands in the central region represent cor- 
related systematic errors related to the uncertainty in 
y m in- As y m in varies the dominant effect is common dis- 
placement of p and q (cf. Fig. HOfl . Multiplicity depends 
mainly on the mode or mean of the fragment distribu- 
tion, and therefore primarily on the ratio or difference 
of p and q. Multiplicity is therefore insensitive to the 
choice of y m in- The FF width on the other hand de- 
pends directly on the sum p + q and is therefore more 
influenced by the choice of y m in- To the right of the left 
dash-dot line the (p, q) vary slowly and linearly with in- 
creasing energy scale. The energy dependence for light 
quarks implies a slight reduction of the mode with the 
peak width unchanged, consistent with the fiducial FFs 
in this study (e.g., Fig. [SJ. The gluon FF shows simi- 
lar mode variation, but the width is also reduced with 
increasing energy. 

Below the left dash-dot line (Q ~ 10 GeV) the (p, q) 
change rapidly. The multiplicity data, especially the 
CLEO data, require a sharp drop in q in that energy 
interval for both quarks and gluons which is effected by 
the tanh term in q q of Eq. JIJ. The convergence of the 
quark and gluon (p, q) at the energy scale defined by the 
lower dotted line, again required by the CLEO data, is ef- 
fected by the tanh terms in p g and q g . Below y m ax = 3.6 
(Q = 5 GeV) there is no guidance from data, but we 
speculate as follows. At 5 GeV the average jet multiplic- 
ity is ~ 2.5 and there is no distinction between quark 



and gluon jets, p ~ q and the FF is therefore symmetric 
about the midpoint on u. We argue that at lower energies 
the mean jet multiplicity approaches one and the FF ap- 
proaches a delta function at u = 1, requiring q — ► 1 and 
p — > oo. We sketch those trends with large error bands 
in the left panel as a simple extrapolation of the trends 
derived from data. 

In Fig. E| (right panel) we show a sequence of model 
functions for nine equal y m ax steps from 1.6 to 7.8, with 
parameters derived from the (p, q) curves in the left panel 
for gluon and quark FFs. The modes for quark jets (solid 
curves) and gluon jets (dashed curves) move from left to 
right with decreasing energy scale, and the dotted curves 
for y m ax < 3.6 (Q < 5 GeV) represent both parton types 
in common. Below 5 GeV the FFs slew to the right and 
may approach a delta- function limit at u = 1. Those 
low-Q 2 trends can be compared with the theoretical de- 
scription of heavy quark fragmentation on the right of 
Fig. El 

The energy dependence of gluon-to-quark-jet multi- 
plicity ratio r = n g /n q derived from the beta parame- 
ters is plotted as the lowest solid curve in Fig. E] (left 
panel) . The open points are taken from [5^] and the solid 
points are from |59j. r is expected to approach the ra- 
tio of color factors Ca/Cf = 2.25 at large Q. The ratio 
indeed increases monotonically with y ma x from unity at 
ymax ~ 3.5 (p ~ 2 GeV/c and n ~ 2), but the approach 
to Ca/Cf is slow. The overall trend is in rough agree- 
ment with theory [571 160| . 

The ratio slope r' = dr/dy max ~ dr/d ln(Q) is also 
of theoretical interest but difficult to calculate since it is 
very sensitive to perturbative corrections |57|, Slope r' 
derived from (p, q) rises to peak value 0.45 at y ma x ~ 
4, then falls to 0.11 at y ma x — 5 (Q = 20 GeV) and 
rises linearly to 0.13 at y max = 8 (Q = 400 GeV). The 
value r' = 0.1 is in rough agreement with theory |5jfl . In 
contrast to the slow evolution of r the rapid separation 
of q q and q g with increasing energy between the lower 
dotted and dash-dot vertical lines of Fig. El (left panel) 
contrasted with a nearly fixed difference between them 
above that region may provide a clearer manifestation of 
the emergence of color charge. 



D. Fragmentation functions on (y,y m ax) 

We can use the parameterized beta distribution 
to construct a 2D fragment distribution on (y, y max ) 
as follows. Form factor (3(u;p,q) describes the 
shapes of FFs over a broad Q 2 interval. The beta 
distribution in turn determines multiplicity n(y max ) 
through (xe) over the same range (Fig. I13fl . We 
combine the two factors to form D(y, y max ) — 
max ) 0[u(y, y m in, Umax); p(y ■max max 1] • 

In Fig. El (left panel) we plot D(y 1 y m ax)- The verti- 
cal dotted and dash-dot lines mark the same energies as 
in Fig. El (left panel). The dashed curve is a 'locus of 
modes' (positions of maxima) of conditional distributions 
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FIG. 15: Left panel: Joint fragment distribution D(y, y m ax) 
on fragment and parton rapidities for inclusive partons (~ 
udsc quarks) and inclusive hadrons. Fragmentation functions 
are vertical slices (conditional distributions) from the joint 
distribution. Right panel: The same distribution transformed 
to (y — y ma x, y max), with y ma x ~ y ~ fp, the logarithmic 
relative momentum. The vertical dash-dot lines define the 
interval determined by fiducial FFs plus dijet multiplicities. 
The intervals between vertical dash-dot and dotted lines are 
defined only by multiplicity trends. The upper-right region of 
the right panel illustrates scaling violations. 



on y for fixed y m ax- The approach of that curve to the 
solid diagonal line (y — y ma x) at lower left corresponds to 
the approach of the dotted curves in Fig. El (right panel) 
to u = 1. The horizontal dotted line denotes y m in, and 
Umax — 8 corresponds to t/s ~ 400 GeV. This joint frag- 
ment density provides the basis for extrapolating FFs 
down to Q ~ 1 GeV (y max ~ 2). Fig. El (right panel) is 
a transformation of the left panel onto (y max , V — Umax), 
with y — y max ~ \nx p = — £ p , which illustrates in the 
upper-right corner scaling violations: variation of the 
fragment density with increasing ln(Q/A) — > y max at 
constant x p or £ p (constant y — y ma x) (cf. Sec. IXIl() . 



XI. PEAK STATISTICS AND PQCD 

We have constructed a simple parameterized model of 
FFs for e + -e~ collisions which compares well with data. 
In this section we compare the energy dependence of peak 
statistics on u and £ p inferred from our parameterization 
with predictions from pQCD. Fragmentation-function 
peak statistics predicted by pQCD [m I27I I6lll6^| can be 
compared to peak statistics u* (mode) , u (mean) and a 2 
(variance) for distribution f3(u;p, q). The mode for the (3 
distribution is u* — p ~ 1 „ and the mean is u — — r- , with 

p+q—2 p+q ' 

(p, q) determined by the parameterizations in Fig. 1141 
The mode on y is y* = u* y max + (1 - u*) y min . 

In Fig. El (left panel) we show measured values of £* in 
the form y max - £* ~ y* (consistent with Fig. gjl vs y max 
[comparable to plots of £* vs ln(Q/A)1 for eig ht quark-jet 
and fourteen gluon-jet energies |2g, [32, |62]- The solid 
curve y*(y m ax) for quark jets inferred from our (p, q) pa- 
rameterization is the same as the dashed curve in Fig. 1151 
(left panel). The five stars are obtained from our fits to 
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FIG. 16: Left panel: Comparison of fragmentation-function 
modes vs parton rapidity from quark and gluon data (points) 
with 'locus of modes' trends (solid and dashed curves) derived 
from (p, q) energy systematics in Fig. 1141 (left panel) and from 
the MLLA (dash-dot and dotted curves). Right panel: Com- 
parison at two energies of FF data, beta distributions on u 
and MLLA gaussians suitably transformed to u. 



the fiducial FFs in Table [U (compare to peak modes in 
Fig. El- The errors are smaller than the points. 

The MLLA predicts for inclusive jets £* = 0.5V 2 + 
C2v^2 - c|, with V 2 = ln(Q/2A) (note the 2 in the 
denominator), c 2 = a/^/16bn c , a — bjr? c and b — 
(lln c - 2n f )/3 31]. The MLLA prediction for £* trans- 
formed to y m ax — £,p is plotted as the dash-dot curve in 
the left panel. The curve corresponds to rif = 5, but 
changes with rif — > 3 are within the data errors. The 
hatched area is the region of interest for study of low- 
Q 2 partons. The MLLA curve diverges from the (p q , g q ) 
parameterization (solid curve) in that region. 

We can also obtain a mode prediction for gluon jets. 
The MLLA prediction for the quark-gluon mode differ- 

ence is A£* = £ - « i (l + + 0(^7) ~ 

0.1 Taking the inclusive £* prediction above as 

£* we plot y max -£,g = y m ax-£,q-A£* as the dotted curve 
in Fig. El which agrees fairly well above y ma x = 4.5 
(Q ~ 12 GeV) with the gluon y* trend (dashed curve) ob- 
tained from parameters (p g , g g ) in Fig.^^ Data from [(j^ 
for FF modes from gluon jets are plotted as solid tri- 
angles. The modes were obtained from gaussian fits to 
gluon FF data over limited intervals on £ p . The data are 
well described by the dashed curve obtained from our 
(p, q) energy systematics and by the MLLA prediction. 

The MLLA width prediction on £ p is a ip = V 3/4 / v / 2cT, 
with ci = ^'iQn c /b |2^|. The width on £ p should be 
equivalent to the width on y (cf. Fig. 0J. The variance 
of the beta distribution on u is err = -, — ; — dr 1 -, — rrr — 

4(p+g+i) ~ 0-035 for flavor-inclusive e + -e~ jets. Thus, 
the observed r.m.s. width on y is a y ~ 0.2y max ~ 0.2V, 
the coefficient nearly independent of y max per the (p, q) 
systematics in Fig. El That result is inconsistent with 
the MLLA width prediction a y - 0.37 V 3 / 4 . 

Measured FFs have been compared directly with ana- 
lytical predictions of peak statistics from the MLLA and 
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with gaussians on £ p denned by parameters from pertur- 
bative approximations ■ I n Fig- El (right panel) we 
compare beta distributions and data for two energies on 
normalized rapidity u with corresponding MLLA gaus- 
sians (normalized to unit integral) using the parameters 
described above. The gaussian tails do not describe the 
data. Our parameterized model is consistent with pQCD 
predictions at larger Q 2 , and the beta distributions (solid 
curves) demonstrate good sensitivity to small but mean- 
ingful systematic variations with energy of the FF data. 
The good fit of beta distributions to data over all frag- 
ment momenta insures a well-defined peak integral. 

XII. SCALING VIOLATIONS 

The naive parton model predicts that parton distribu- 
tion functions (PDFs) measured in deep inelastic scatter- 
ing [analogous to FFs D(x, Q 2 )} should be independent of 
energy scale Q 2 (Feynman-Bjorken scaling) |6ij. Scaling 
violations |yj — variations of PDFs and FFs with energy 
scale — are described by the DGLAP equations [33, |33 
and depend on the running of strong coupling constant 
a s , on the available phase space for fragmentation and on 
1 — > 2 parton splitting described by the Altarelli-Parisi 
splitting functions. Scaling violations of measured FFs 
can in turn be used to determine a s (cf. Fig. 118(1 and to 
test the predicted values of QCD color factors Ca and 
Cp. In this section we describe scaling violations on con- 
ventional momentum/energy fractions and on rapidity 
variables in terms of the FF parameterization developed 
in this study. 

A. Scaling violations on (x,Q 2 ) 

Fragment distributions on x p are approximately expo- 
nential (cf. Fig. [21 - left panel). Scaling violations are 
described in that format as follows. The slope of the FF 
on x p becomes more negative (the distribution 'softens') 
with increasing energy scale. Scaling violations (slope 
change with increasing energy scale) are larger for gluon 
jets than for quark jets (more radiation is produced by 
the larger effective color charge of gluons) . The strength 
of scaling violations at large x p and the accompanying 
increase of the FF at small x p (another manifestation of 
scaling violations) are directly related by energy conser- 
vation, as noted in Sec. IIXI Comparisons of quark-jet 
FFs at different energy scales have been used to measure 
the running of a s [28l Eol l64j . Comparisons of gluon- 
and quark-jet FFs have been used to measure the ratio 
of color factors C A /C F jH |H |H En- 
sealing violations are described by QCD theory in the 
form of the DGLAP equations 38, 39], which are in lead- 
ing order (LO) 

= ^ rU_z 

dins 2n ' J x z 



P a b(z) are the Altarelli-Parisi splitting functions |3^|, and 
a, b denote parton combinations. In a typical study of 
scaling violations FFs are parameterized at several en- 
ergy scales s with a model function such as D(x, s) — 
Nx a (1 - xf (1 + 7/ar) [5J IU Such parame- 

tcrizations can be quite extensive. The KKP parame- 
terization j4lj | employs 14 parameters for each parton- 
hadron combination, the energy dependence of each of 
(iV, ct,/3, 7) being described by several polynomial coef- 
ficients. The parameters are determined by using the 
DGLAP equations to evolve the model FFs across en- 
ergy scales, varying (N, a, (3, 7) with energy to best fit the 
data and emphasizing the region x > 0.1 where pQCD is 
most applicable (cf. Fig. El - left panel). 

To illustrate scaling violations in a conventional 
context with the results of the present study we 
transform parameterized joint fragment distribution 
D(y,y max ) in Sec. KEI (Fig- 113 to D(x E ,Q 2 ) = 
p/(Ex E )D[y(x E ,Q),y max (Q)]. In Fig. El (left panel) we 
plot conditional distributions D(x E ,Q 2 ) for x E = 0.02, 
0.07, 0.15, 0.27, 0.41, 0.60, 0.81 vs Q = m cosh(y max ). 
The curves for both udsc jets (solid) and gluon jets 
(dashed) compare well with previous analyses (e.g., Fig. 
10 of 46]). The general trends agree with the descrip- 
tion of scaling violations noted above but extend over a 
broader energy range than is usually obtained from data. 
The sharp falloffs at smaller Q and x E occur at kinematic 
limits lax E ~ y m m — ymax defined by the dotted line in 
Fig. El (right panel) . Other features of the distributions 
correspond to structures in the (p, q) trends of Fig. 1141 
(left panel). In Fig. 1171 fright panel) we replot the same 
curves vs y m ax on a logarithmic scale. Above y max = 5 
(Q = 20 GeV) the curves are nearly straight, revealing 
the power-law behavior expected for pQCD. 




1 10 10 2 2 

Q (GeV) y, 



FIG. 17: Scaling violations in a conventional format. Each 
curve is a conditional slice at fixed xb from the distribution 
2dn/ dxE(xE,Q 2 ) suitably transformed from model function 
2n(y max ) fi(u(y);p,q) in Fig. E] (right panel). 

The vertical dotted lines in both panels separate three 
regions. Region A (Q = 1 — 5 GeV) is dominated by non- 
perturbative effects. Color charge is effectively hidden 
(quarks appear similar to gluons) and parton fragment 
multiplicities are ~ 1 - 3. Although it is least amenable 
to theoretical treatment, region A produces the major- 
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ity of partem fragments in nuclear collisions and there- 
fore requires at least a phenomcnological characteriza- 
tion consistent with QCD theory. Extrapolation of FF 
systematics into this region is the purpose of the present 
analysis. 

Region B (Q = 5 — 20 GeV) is the transition region in 
which color emerges and fragmentation approaches the 
perturbative description. Convergence of the gluon FF 
with the quark FF near and below Q = 10 GeV in Fig. 1171 
determined in this analysis by multiplicity trends alone, 
is also apparent in direct measurements of scaling vio- 
lations in FFs 46], for example, deviation of the gluon 
FF from the HERWIG Monte Carlo at large x E and for 
energy scale (Q jc t = Q/2) 5 - 10 GeV in Fig. 6 of 46]. 

Parameters (p, q) vary weakly and linearly in the en- 
ergy range above Q — 20 GeV (cf. Fig. E]~ left panel). 
If we set the slopes of p and q on y ma x to zero in that 
interval the changes in Fig. El are small compared to 
the dominant structure. We conclude that much of the 
variation in the perturbative region of Fig. El is deter- 
mined by phase-space acceptance variations with parton 
energy. The subtle linear variations of the (p, q) param- 
eters in that region may provide more differential access 
to the parton cascade process. To explore that possibility 
we consider a modified form of the DGLAP equations in 
the remainder of this section. 



B. The running of a s (Q) 

The energy-dependent a s factor in the LO DGLAP 
equations can be approximately eliminated as follows 
(this procedure reverses use of the energy dependence 
of scaling violations to infer a s 28]). Data showing the 
running of a s with en ergy scale Q are summarized in 
Fig. El (left panel) |65t log using a conventional plotting 
format. A NLO expression for a s (Q) is 67] 



a s (Q) 



2tt 



1 



fa ln(2V) 
0{ 2Y 



■0(1) 



{2Yf 



(3) 



with Y = ln(Q/A), fa = (lln c — 2n/)/3, n c and n/ the 
color and effective flavor numbers and fa — (102n c — 
38n/)/3. Eq. is plotted as the dash-dot curves in 
Fig. ED For A = 0.2 GeV, n f = 5 and 0(1) = 0.3 the 
NLO curves describe the data well. The solid curve in 
the left panel is obtained from the straight-line parame- 
terization in the right panel. 

Fig. El ( r ight panel) shows an alternative plotting for- 
mat in which the data are consistent with a linear trend 
on rapidity, and the low-Q 2 region is visually more acces- 
sible. In the leading-log (LL) approximation the strong 
coupling constant is given by |39j 



l/a s (Q) = l/a s (Q ) + fa/2ir ■ \n(Q/Q ) 



(4) 



Since y m ax{Q A) ~ ln(Q/A) for Q>Awe parameter- 
ize a s in the form l/a s (Q) — A(A) y max (Q/2; A), with 
the constraint that for any A the line should pass through 
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FIG. 18: Left panel: Conventional presentation of a s sys- 
tematics on energy scale Q. Right panel: Systematics of l/a s 
on rapidity y ma x- The dotted lines cross at a 3 (Mz )- 



the current average a s (Mz ) — 0.118 (crossed dotted 
lines). The linear relation is then a single-parameter 
expression since A(A) w [l/a a (M Zo )}/ln(M Zo /A). The 
range of A values permitted by one-sigma deviations at 
both the 1.78 and 4.1 GeV data points is A = 0.2 ± 0.05 
GeV, consistent with the adopted NLO value above. The 
corresponding A value is A — 1.37 ± 0.05, which can be 
compared with the LL slope fa/2ir = 1.22 for n/ = 5 
flavors. The hatched box denotes the region of primary 
interest for the study of low-Q 2 parton fragmentation 
(ce s 6 [0.2,0.4]). The approximate proportionality of 
l/a s and y m ax apparent in both the LL and NLO de- 
scriptions is used in what follows. 



C. Scaling violations on (y,y m ax) 

We have determined that D(x,s) — ► D(£, Ins) — > 
D(y,y max ) has a simple underlying structure (beta dis- 
tribution) and energy dependence [{p, q) parameters]. 
What are the implications for pQCD and the DGLAP 
equations? The general behavior of joint distribution 
D(y,y m ax) is evident from the surface plots in Fig. 1151 
where the trivial £ term in hiD(x.s) = ]nD(£, s) + £ 
which dominates Fig. 1171 is eliminated, permitting more 
precise comparisons on a linear scale. Eq. (J2J has the 
form dD(x, s)/dhis oc a s (s) x convolution integral. The 
running coupling constant introduces an energy depen- 
dence which can be eliminated. 

As an alternative approach we introduce the logarith- 
mic derivative |62| motivated by the relation between 
Mcllin transforms of FFs and splitting functions. The 
Mcllin transform of an FF is 



D{w,s) 



dxx™- 1 D(x,s) 



(5) 



which is also the Laplace transform of D(^,s). The 
DGLAP equations, written in terms of Mellin transforms, 
are represented by a simple matrix equation |39|] . For the 
non-singlet case D ns = D„ — Dq the logarithmic deriva- 
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tive is 



dhiD DS (w,s) a s (s) 



d In 



2tt 



Pq q (w)=J qq (w,s), (6) 



where P qq (w) is the Mellin transform of splitting func- 
tion P qq (z) for the process q — * q(z) + g(l-z), with 
the respective momentum fractions noted, and "f a b(w,s) 
are the anomalous dimensions of QCD |68l l69| . Since 
dins « 2dy max we multiply through by 2y max and use 
the results of the previous subsection to obtain 



dlnD ns (w,y max ) 1 
d\ny max ttA 



(7) 



For the non-singlet case the logarithmic derivative of the 
Mellin transform of an FF is proportional to the Mellin 
transform of a splitting function, independent of energy 
scale in LO. That simple relation motivates a similar ap- 
proach to the FFs themselves. 

We multiply Eq. J2J) by xhx{s/rno)/xD{x, s) to obtain 
dlnD(f,s)/dln[ln(s/mg)] -> dhxD{y, y max )/dln 
on the LHS. The additional factor ln(s/mg) ~ 2y max 
introduced in converting to the logarithmic derivative 
cancels the l/y m ax trend of a s (y ma x) on the RHS to 
good approximation (Sec. IXII B|) . What remains on 
the RHS is the convolution integral, including split- 
ting function zP(z) = P(Q with £ = ln(l/z) and 
fragmentation-function ratio D a (£ — £, s)/JD(,(£, s). Wc 
then use the following transformations, £ — > y max — y, 
C ~> Umax - y', D(£) -> D(y max - ~ D (u) and 

D(£ - C) -» D(y max - £ + ~ D (V + Umax ~ y'), to 

obtain 



din D b (y,ymax) 

d 111 Umax 



a 

D a (y + y, 



dy ' P a b(y m ax - y') (8) 



Db(y, ymax) 

a modified form of the DGLAP equations on (y, y m ax)- 

In Fig.HH(left panel) we plot d\nD(y,y max )/d\ny max 
vs y m ax for quark (solid) and gluon (dashed) FFs us- 
ing the parameterized D(y,y max ) from the present anal- 
ysis for each parton type. The structure is completely 
defined by parameters (p, q) in Fig. [21 (left panel) and 
corresponds exactly to the scaling violations in Fig. El 
Because of the form of the beta distribution the loga- 
rithmic derivative directly relates exponents p and q to 
the splitting functions. There are three main features 
in the distributions: 1) nearly-horizontal linear trends at 
larger energy scales (to the right of the dotted line), 2) 
'singularities' at smaller energy scales due to kinematic 
boundaries and 3) minima at intermediate energies cor- 
responding to a transition from 'small' (1-2) to 'large' 
(3 or more) fragment number, which may also relate to 
the emergence of color charge (quark-gluon distinction) 
at Q ~ 8 GeV. The nearly constant values for larger 
xe (y —> ymax) and y max are the slopes or 'power-law' 
exponents in Fig. 1171 fright panel). 
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FIG. 19: Left panel: Scaling violations in the form of loga- 
rithmic derivatives of the distributions in Fig. 1171 (left panel) . 
Solid curves represent udsc quark jets, dashed curves repre- 
sent gluon jets. The near uniformity to the right of the dotted 
line (Q — 20 GeV) for larger xe is notable. Right panel: Log- 
arithmic derivatives for quark and gluon jets at y m ax = 7.8 
vs energy fraction xe = cosh(y)/ cosh(y ma x)- 



In Fig. UU (right panel) din D{y,y max )/d\n 

ymax VS 

xe = cosh(y)/ cosh(y max ) is plotted for quarks and glu- 
ons. In the limit xe — ► 1 (y — > ymax) and large y, ™° 
expect |62] | 



we 



{d\nD(y,y max )/d\n ymax} gluon C A 
{d In D(y, y ma x)/d In y max } 

quark 



= 2.25. 



(9) 



The dotted lines are o(xe + b) and 2.25 x o(xe + b), 
with a = —1.8 and b = 0.35 adjusted to best match the 
quark points. The ratio trend is in reasonable agreement 
with the QCD expectation above xe — 0.2. Within the 
constraints established in this analysis there is freedom 
to adjust the energy dependence of (p g ,q g ) to achieve 
that agreement. Changing the slopes of (p s , q g ) on y max 
changes the position of the gluon curve in Fig. ^5] (right 
panel). Variations of ~ 20% sketched by the lower 
hatched area are possible without disturbing agreement 
with data. 

The opposite-sign slopes for (p q , g q ) in Eq. Q are 
strongly constrained by the quark-jet multiplicity trend 
2n q (y max ) in Fig. ^| Setting those slopes to zero gives 
the lower dotted (quadratic) curve in the right panel of 
that figure expected for self-similar FF scaling and fixed 
mode u*. Varying the same-sign slopes in the (p g ,g g ) 
expressions together does not affect the gluon-jet multi- 
plicity trend 2n g (y max ) determined by the q— p difference 
but does affect the slopes of the gluon logarithmic deriva- 
tives at larger xe in Fig. ^5] (left panel) and the slope of 
the gluon curve in the right panel. Setting the same-sign 
(Pgj Qg) slopes to zero makes the logarithmic derivative 
slopes zero, and the gluon curve in the right panel be- 
comes parallel to the quark curve (upper edge of lower 
hatched region) . Increasing the (p g ,q g ) same-sign slopes, 
which reduces the width of the gluon FF, and hence its 
amplitude at larger xe, without changing its mode or 
the multiplicity trend provides a match to the expected 
Ca/Cf ratio. 
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XIII. DISCUSSION 

Our intention in this study has been to provide the phc- 
nomenological means to extrapolate e + -e~ parton frag- 
mentation functions to low Q 2 where the perturbative de- 
scription of QCD is not applicable. We have been guided 
by the many precise measurements of e + -e _ fragmen- 
tation functions and their relationship to pQCD predic- 
tions now available. Parton scattering and fragmentation 
at low Q 2 are in turn important for understanding p-p 
and A- A collisions at RHIC. As a result of this study we 
have found that the beta distribution provides a simple 
but precise description of FFs which accomplishes the de- 
sired extrapolation but also reveals some interesting new 
aspects of parton fragmentation. 

A. Motivation from p-p and A-A collisions 

We observe large-amplitude two-particle correlations 
in nuclear collisions at RHIC, produced in part by frag- 
mentation of low-Q 2 partons (minijets). Minijet-related 
correlations observed in p-p collisions are strongly modi- 
fied with increasing centrality in Au-Au collisions. While 
low-Q 2 partons play an important role in forming the col- 
ored medium and driving large-scale hydrodynamic phe- 
nomena according to theory, they may also function as 
sensitive probes of that medium. However, theoretical 
descriptions of low-Q 2 scattering and fragmentation are 
limited. Factorization is not applicable since low-Q 2 par- 
ton scattering and fragmentation remain intimately con- 
nected. However, new aspects of fragmentation observed 
via correlations in nuclear collisions (including strong de- 
pendence of fragment angular correlations on Q 2 0) sug- 
gest a complex but understandable low-Q 2 process. 

To facilitate theoretical descriptions of low-Q 2 phe- 
nomena we have attempted to extrapolate a phenomeno- 
logical representation of measured FFs in e + -e~ collisions 
to low Q 2 . We want to connect two-particle fragment 
correlations in p-p and heavy ion collisions to pQCD and 
conventional jet phenomenology through single-particle 
fragmentation functions. The extrapolation imposes spe- 
cial demands on the fragment representation (particu- 
larly for small particle momenta) which have led us to 
employ rapidity y and normalized rapidity u as our basic 
kinematic variables. That decision led to the discovery 
that the beta distribution on u is a good model of light- 
quark and gluon fragmentation functions. 

In p-p and A-A collisions at RHIC we encounter co- 
pious parton fragmentation in an energy regime where 
pQCD is not applicable, but where trends from pQCD 
may provide semi-quantitative guidance for analysis and 
interpretation. We therefore distinguish between condi- 
tional pQCD fragmentation functions and unconditional 
fragment distributions measured in nuclear collisions. 
Given that distinction we can attempt to connect low- 
Q 2 phenomena in nuclear collisions to QCD through the 
close connection between FFs and fragment distributions 



as a limiting case. This is a new aspect of fragmentation 
which lies outside the scope of conventional pQCD frag- 
mentation analysis. 

B. The beta distribution as compact representation 

FFs for e + -e~ collisions transformed to y are approx- 
imately self-similar with increasing y m ax, as sketched in 
Fig. n They are bounded by parton rapidity y max and 
lower limit y m in- Those trends suggest a further trans- 
formation to normalized rapidity u — (y — ymin)/ (y max — 
y m in)- Measured FFs plotted on u are nearly indepen- 
dent of parton energy and can be factorized into di- 
jet multiplicity 2n(y max ) and unit-normal form factor 
g(u, y m ax) which is modeled by the beta distribution. 
Since energy conservation relates the dijet multiplicity 
to the form-factor shape the fragmentation process is 
completely represented by the energy-scale dependence 
of parameters (p, q) of the beta distribution. 

The beta distribution with two energy-dependent pa- 
rameters thus precisely describes light-flavor FF data 
over the scale interval Q G [5, 100] GeV (e.g., Figs.EJ-EJ) 
and can be extended to lower energies, accomplishing the 
main goal of this study: extrapolation of fragmentation 
systematics down to fragment multiplicity n ~ 1-2 and 
energy scale Q ~ 1 GeV. We also obtain simple repre- 
sentations of FFs on (x,s) or (y,y m ax) (continuous 2D 
surface) for studies of scaling violations over the full kine- 
matic range. 

C. Fragmentation as an equilibration process 

Why does the beta distribution provide a good de- 
scription of e + -e~ FFs? (3(u;p,q) describes systems in 
which entropy is maximized (e.g., by a parton cascade) 
on a bounded interval (e.g., bounded by the leading- 
parton energy). /3(u;p,q) is one instance of the ex- 
ponential family of probability distributions p(x) which 
can be defined in terms of a maximum-entropy condi- 
tion with constraints [jjj- The beta distribution de- 
fined on a; G [0, 1] maximizes the Shannon entropy 
S = — J dxp(x) \n[p(x)\ subject to constraints on geo- 
metric means ln(ai) = J dxp(x) \n(x) and ln(l — x). The 
gaussian distribution can be similarly defined, with con- 
straints on its first and second moments x and x 2 . 

Given those properties of the beta distribution and its 
good description of FFs, fragmentation of light quarks 
and gluons can be viewed as an equilibration process 
controlled by two opposing tendencies: parton splitting 
as a form of downscale energy transport which increases 
entropy and gluon coherence which constrains the split- 
ting at a scale conjugate to hadron size. The observed 
fragment distribution is then a maximum-entropy con- 
figuration balancing those two tendencies. Of the two 
beta parameters q reflects the splitting tendency and p re- 
flects the hadron size constraint. The DGLAP equations 
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which describe the perturbative splitting (transport) pro- 
cess are thus coupled to the soft part of the FF by entropy 
maximization as well as energy conservation. 



a context for the present phenomenological analysis. 



E. Energy dependence and scaling violations 



D. Conventional pQCD fragmentation functions 

Do benefits from the introduction of (y,y m ax), u and 
the beta distribution justify a new approach to fragmen- 
tation? In Fig. 1201 we compare FF data to a beta distribu- 
tion fit from this analysis and a pQCD model FF (KKP) 
obtained from a conventional scaling violations analysis 
using the DGLAP equations ^lj an( i defined by 14 pa- 
rameters for each parton-hadron combination. Fig. 1201 
(left panel) shows the OPAL 91 GeV FF data from Fig.H 
with the KKP fragmentation function (dashed curve) and 
the FF from this analysis (solid curve) . With the excep- 
tion of a small deviation at large x p the agreement in this 
format appears to be good (also cf. Fig. 0- left panel). 
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FIG. 20: Left panel: Beta distribution (solid) and KKP FF 
(dashed) curves compared to OPAL 91 GeV data points (open 
circles) on linear momentum variable x p . Right panel: The 
same curves and data transformed to normalized rapidity u. 
The vertical dotted lines both correspond to x p = 0.1. 

Fig. EDI (right panel) shows the same distributions on 
normalized rapidity u. The KKP FF based on DGLAP 
evolution deviates strongly from data below u = 0.7 and 
accurately represents less than 10% of the fragments at 
91 GeV. The vertical dotted line in each panel shows 
the intended region of validity (x p > 0.1) of the KKP 
and similar FFs. In contrast, the FF from our analysis, 
based on the beta distribution and defined by two param- 
eters per parton-hadron combination with simple energy 
dependence above Q — 20 GeV, accurately describes the 
data on x p over six orders of magnitude and, extrapolat- 
ing the full data distribution down to zero momentum, 
provides a well-defined multiplicity integral. 

Recent improvements in the pQCD description of FFs 
have lead to much improved coverage on x p . In |24| a 
consistent combination of DGLAP evolution, resumma- 
tion of soft gluon logarithms and incorporation of hadron 
mass effects provides a semi-quantitative description be- 
low the FF mode on y or u while retaining good agree- 
ment above the mode. Such theoretical advances provide 



We have studied the energy dependence of FFs as rep- 
resented by beta distribution parameters (p, q) in Figs. 1 131 
1161 The direct relation between dijet multiplicities and 
FF shape parameters via the energy sum rule provides 
a new method for extending the FF description over a 
broad energy range. The connection between shape pa- 
rameters (p, q) and dijet multiplicity 2n(y max ) is partic- 
ularly important for extrapolation to low Q 2 (e.g., the 
CLEO multiplicity data provide strong constraints on 
low-energy FF evolution) . 

A striking feature of Fig. ^|is the bifurcation of quark 
and gluon trends near Q ~ 5 GeV and evolution to large 
separation by Q ~ 20 GeV, with a dramatic correspon- 
dence in Fig. E| (left panel). Above 20 GeV the system 
exhibits simple pQCD trends which are most apparent 
in Fig. 1141 Corresponding peak modes on u remain close 
to 0.4 over a broad energy range. It is not clear from 
Fig.^|that peak modes increase rapidly near 5 GeV and 
approach unity below that point, but that trend can be 
inferred from Fig. 1141 which provides important extrapo- 
lation guidance. Mode trends are illustrated in Fig. ED 
and agree well with data and pQCD predictions in the 
energy range above 15 GeV. 

We have considered scaling violations in Fig. ED and 
Figs. We find that a combination of QCD scaling 
violations and gluon coherence leads to approximate self- 
similarity of FFs on rapidity with changing energy scale, 
as illustrated in Fig. E and demonstrated in Fig. [3J Nor- 
malized FFs plotted on normalized rapidity u are there- 
fore nearly independent of energy scale. However, rela- 
tive to the self-similarity trend we observe more subtle 
forms of 'scaling violations' as described in the previous 
paragraph. 

Conventional scaling-violation systematics are easily 
and precisely reproduced by our parameterization over a 
broad energy range, as demonstrated in Fig. 1171 down to 
kinematic limits. It is straightforward to explore the con- 
sequences of varying (p, q) energy trends. For example, 
Fig. 1191 fright panel) demonstrates that the Ca/Cf limit 
for logarithmic derivatives previously establishe d by spe- 
cific experimental measurements (e.g., pa. I54L l59l \62\ ) 
is consistent with the (p, q) parameterization determined 
by the present study. 



XIV. SUMMARY 

In conclusion, we have examined the variation of frag- 
mentation functions (FF) with parton energy on conven- 
tional kinematic variables x p and £ p , on rapidity y and on 
normalized rapidity u. We find that FFs plotted on ra- 
pidity y vary with energy in a nearly self-similar manner. 
FFs transformed to u are well described by a product of 
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the dijet multiplicity and a unit-normal form factor mod- 
eled by the beta distribution. The latter is determined by 
parameters (p, q) which exhibit modest linear variations 
within perturbative energy scale range Q > 20 GeV. The 
beta distribution shape, when combined with an energy- 
conservation sum rule, also determines FF multiplicities. 
The factored representation on u thus provides a simple 
and compact representation of e + -e~ FFs over a broad 
energy range and permits extrapolation to small energy 
scales. 

Reduction of scaling violations to near-linear variations 
of beta parameters (j>, q) may provide more differential 
access to the energy dependence of parton fragmenta- 



tion. The beta distribution model suggests that frag- 
mentation of light partons can be viewed at larger energy 
scales as an entropy-maximizing equilibration process. In 
this analysis the energy dependence of fragmentation has 
been extrapolated down to a kinematic region not ac- 
cessed by conventional methods. Such low-Q 2 extrapo- 
lation provides a phcnomcnological context for minijct- 
related two-particle correlations in p-p and A-A collisions 
at RHIC, forming a basis for theoretical treatments of 
in-medium dissipation of low-Q 2 partons and the subse- 
quent hadronization process in heavy ion collisions. 

This work was supported in part by the Office of Sci- 
ence of the U.S. DoE under grant DE-FG03-97ER41020. 
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